(  10  ) 


where  fl  Is  the  rotor  RPM  change  from  trim  sensed  by  the  governor 
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HOVERING  HELICOPTER  FLIGHT  DYNAMICS  - 
A  STUDY  OF  VERTICAL  MOTION  (U) 


by 

C.R.  GUY 


SUMMARY 

I 

This  document  examines  the  flight  dynamic  characteristics  of 
a  hovering  helicopter  in  response  to  collective  stick  inputs.  Only 
vertical  motion  is  considered.  A  mathematical  model  of  the 
aerodynamics/kinematics  and  rotor  system  dynamics  of  a 
representative  aircraft,  linearized  about  the  hovering  flight 
condition,  is  used  as  a  basis  for  the  study.  The  dynamic  behaviour 
of  this  model  is  compared  with  flight  test  results  for  similar 
inputs  and  the  coefficients  are  adjusted  to  improve  model 
performance.  The  model  structure  is  then  modified  by  the 
incorporation  of  additional  terms  to  improve  its  accuracy.  Although 
no  attempt  is  made  to  optimize  the  coefficients  of  the  model  using 
parameter  identification  techniques,  the  approach  used  here  provides 
acceptable  results.  - 
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1 .  INTRODUCTION 


(  1  ) 


The  complex  aerodynamics  and  flight  mechanics  associated  with 
helicopters  makes  accurate  mathematical  modelling  of  flight  dynamic 
character lsti cs  difficult.  The  behaviour  of  such  models  when  compared 
with  flight  test  results  for  similar  inputs  often  only  indicates  the 
general  trend  of  response.  Some  examples  of  such  models  and  comparisons 
of  their  behaviour  with  flight  test  results  are  given  in  references  1  to 
4.  The  models  used  are  usually  relatively  large,  non-linear  and 
complicated,  thus  making  them  unsuited  to  techniques  such  as  parameter 
identif lcation  which  might  be  used  to  obtain  improved  representation. 

In  order  to  understand  helicopter  flight  dynamic  characteristics 
more  precisely  and  obtain  improved  mathematical  models,  it  is  necessary 
to  simplify  the  task  by  subdividing  it.  As  the  hovering  flight  case  is 
perhaps  the  most  straightforward  to  consider,  only  this  part  of  the 
flight  envelope  is  addressed  here.  In  addition,  this  document  considers 
only  vertical  motion  response  to  collective  stick  Inputs  without 
interactions  in  other  motions.  However,  similar  techniques  to  those 
developed  here  can  be  used  to  determine  the  model  structure  for  yawing, 
pitching  and  rolling  motions  as  well  as  interactions,  as  demonstrated  in 
reference  5.  Flight  test  results  illustrating  dynamic  response  to  stick 
and  pedal  inputs  for  a  Sea  King  helicopter  in  this  condition  are 
available  (reference  6),  as  is  a  flight  dynamic  mathematical  model 
linearized  about  the  hover  condition  (reference  7).  Although  the  flight 
test  results  and  model  relate  to  different  versions  of  this  aircraft, 
the  basic  model  structure  should  he  applicable  to  both  versions. 


A  linear  model  rather  than  a  more  complex  non-1 Inear  one  (e.g.  as 
described  In  reference  4)  was  selected  both  to  simplify  the  problem  and 
understand  Its  fundamental  nature.  In  addition,  a  linear  model  lends 
itself  to  the  application  of  parameter  Identification  techniques  at  a 
later  stage  to  obtain  optimum  fitting  of  the  model  to  flight  data.  The 
linear  hover  model  would  be  obtained  from  the  six  degree  of  freedom  non¬ 
linear  equations  of  motion  using  small  perturbation  theory.  As  the 
manoeuvres  considered  are  essentially  small  perturbations  from  the 
steady  state  hover  condition,  the  linearized  model  approach  is 
justified. 

The  exercise  is  essentially  one  entailing  fitting  of  coefficients 
within  an  existing  model  structure  followed  by  modification  of  structure 
through  the  addition  of  terms.  A  comparison  of  model  behaviour  with 
flight  test  results  for  the  3ame  input  is  made  at  each  3tage  to  ensure 
reasonable  agreement  is  achieved.  In  adopting  this  approach  it  was 
found  necessary  to  add  terms  to  the  model  which  were  not  included  in  the 
basic  model  of  reference  7. 

In  determining  vertical  motion  of  the  aircraft  about  the  hover,  two 
separate  'sub-systems'  influence  overall  behaviour.  These  are  the  rotor 
system  dynamics  and  the  aircraft  dynamics.  The  mathematical  modelling 
of  the  various  elements  of  each  'sub-system'  is  described  together  with 
simulation  and  some  analytical  studies.  Comparisons  of  model  behaviour 
with  flight  test  results  are  also  presented,  along  with  some  notes  on 
stabilizer  systems. 


2 


ROTOR  SYSTEM  DYNAMICS 


2.1.  General 

Flight  teat  results  from  a  Sea  King  helicopter  (reference  6) 
indicate  that  variations  occur  in  engine  torque  and  rotor  speed  when  a 
collective  stick  input  is  applied  to  the  aircraft  in  steady  hovering 
flight  (e.g.  Figure  1).  It  also  appears  that  the  normal  acceleration  of 
the  aircraft  contains  an  oscillatory  component  related  to  engine  torque 
variations.  For  these  reasons,  and  the  need  to  consider  torque 
variations  when  analysing  yaw  response,  the  rotor  system  dynamics  need 
to  be  included  in  the  flight  dynamic  mathematical  model  of  the 
aircraft.  Earlier  helicopter  mathematical  models  (e.g.  reference  8) 
assume  constant  rotor  speed  during  manoeuvres  whereas  more  recent 
studies  (e.g.  references  3  and  4)  have  taken  rotor  system  dynamics  into 
account. 

The  main  purpose  of  this  work  is  to  determine  a  realistic  model 
structure  that  gives  representative  dynamic  behaviour  of  the  system 
being  considered.  This  can  then  be  used  for  stability  analysis  and 
simulation  purposes.  In  line  with  the  configuration  of  the  aircraft 
vertical  motion  dynamics  (section  3)  which  is  based  on  reference  7,  a 
linear  model  of  the  form: 

X  ■  Ax  +  By  ( 1  ) 

is  used  for  the  rotor  system  where 

x  is  the  state  vector 


ft 


y  Is  the  control  vector 


A  is  the  state  matrix  and 

B  is  the  control  matrix. 

As  stated  previously,  the  approach  used  derives  the  model  and  its 
coefficients  from  an  amalgamation  of  established  theory  and  flight  test 
results . 

Instrumented  variables  available  from  the  flight  test  results  of 
reference  6  which  are  involved  in  the  rotor  system  dynamics  mode]  are 
collective  stick  position,  collective  blade  angle,  engine  torque  and 
rotor  speed  (RPM).  Parameters  have  been  determined  so  that  a  reasonable 
approximation  to  dynamic  response  is  obtained,  but  no  attempt  has  been 
made  at  this  stage  to  optimize  the  coefficients. 

2.2.  Mathematical  Modelling 

Based  on  information  from  a  number  of  sources,  a  block  diagram  of 
the  elements  forming  the  simplified  rotor  system  dynamics  for  a  Sea  King 
helicopter  can  be  constructed  as  shown  in  Figure  2.  A  breakdown  of  the 
elements  shown  on  this  diagram.  Including  mathematical  representations, 
is  given  below. 

It  should  be  noted  that  the  model  developed  here  is  linearized 
about  a  specific  flight  condition.  In  this  case  hovering  flight  at  a 
given  aircraft  weight.  As  shown,  transfer  functions  and  equations 


relate  to  changes  in  variables  about  this  flight  condition.  Hence  to 
get  actual  values,  constants  and  integrator  initial  conditions  need  to 
be  added  where  appropriate.  These  are  considered  in  Section  2.3- 


2.2.1  ,  Flying  controls 

A  detailed,  non-linear  mathematical  model  of  the  flying  controls 
for  the  Sea  King  helicopter  has  been  described  in  reference  9.  However, 
for  tnis  purpose  it  is  sufficient  to  use  a  linear  representation  of  the 
form  adopted  by  Biakeiocn  (reference  10),  witn  time  constant  and  gain 
determined  through  ground  tests  on  a  Sea  King  aircraft.  In  transfer 
function  form  this  is: 


K, 

1  +  T  ,  3 


(2) 


or,  in  equation  form  (omitting  initial  conditions): 


K,  , 

®C  ■  t  77  ^  0CST  t  —  J  0C 


where  8C  «  collective  blade  angle  ;  change  from  trim 
9^  -  collective  stick  angle;  change  from  trim 

K1  =  collective  flying  controls  gearing 

t,  -  collective  flying  controls  time  constant 


s 


=  Laplace  operator  and 


2.2.2.  Collective  pitch  anticipator 


Reference  11  states  that  the  output  signal  from  the  collective 
pitch  anticipator  comprises  two  terms  thus: 

i.  A  rate  component  which  is  used  as  feedforward  to  enhance  the  engine 
response  to  changes  in  collective  pitch:  i.e.  it  helps  reduce 
transient  droop,  and 


ii.  A  position  component  which  moves  the  droop  law  thereby  producing 
overcancelled  static  droop. 

These  characteristics  can  be  modelled  by  a  transfer  function  of  the 

form: 


nREF  Kj(l  +  TzS) 

®CST  ’  (1  *  T’s) 


K2 

—  )  e„ 


®CST  ^  t,  ^  nREF 


(«) 

(5) 


where  SI,,-,,  =  rotor  reference  RPM;  change  from  trim 

n  Lr 

Kz  *  anticipator  gain  and 
t2,t3  *  anticipator  time  constants 


Values  for  K2 ,  t2  and  t,  were  found  from  the  simulation  by  trial 
and  error,  flight  trial  and  model  results  being  compared. 


2.2.3.  Governor 


The  governor  compares  the  rotor  RPM  relative  to  the  fuselage  (ft) 
with  the  rotor  reference  RPM  (ftD_„)  to  produce  the  RPM  error 

ntf 

signal  (Rg).  The  demand  signal  to  the  fuel  flow  computer  is 
proportional  to  ftg.  Consequently: 

qe  ■  qref  '  fl  C6) 

2.2.4.  Fuel  flow  coaputer 

The  fuel  flow  computer  produces  a  fuel  flow  change  related  to  the 
demand  (RPM  error)  signal  ft^,.  Padfield  (reference  3)  has  indicated  that 
the  transfer  function  for  this  is: 

«f  k3 

ftg  "  1  ♦  T„S 

K,  , 

i.e.  «f  -  [-)  «E  -  (  —  )  «f 

where  Uj.  »  fuel  flow  rate;  change  from  trim 
K,  >  computer  gain  and 
t,  *  computer  time  constant 

The  gain  and  time  constant  were  determined  from  a  comparison  of 
model  run3  with  flight  test  results. 


(7) 

(8) 


2.2.5.  Qigine 


Reference  3  also  describes  the  relationship  between  engine  torque 
change  from  trim  (Qg)  and  fuel  flow  rate  change  (u^.)  to  be: 

Q  K„ ( 1  +  t,s) 

=;  ‘  -  (i  *  TTsT  "  (9) 

i.e.  Qe  -  (  ^  )  u,  ♦  (  )  ;  -  (  -1-  )  Q  (10) 

t  T,  I  Tt  I  T6t 

where  K„ ,  ts  and  ts  are  gains  and  time  constants.  The  reference  states 
that  r5  and  t6  vary  with  QE  (or,  more  correctly,  high  speed  turbine  RPM) 
as  illustrated  in  Figure  3.  However,  from  the  flight  results  of 
reference  6,  the  overall  system  response  did  not  appear  sensitive  to 
variations  in  t5  and  t*.  Consequently,  in  order  to  retain 

linearity,  t5  and  t6  were  chosen  to  be  constant  and  were  net  to 
correspond  approximately  to  the  hover  torque  value  at  a  specified 
aircraft  weight. 

2.2.6.  Rotor 

The  rotor  torque  change  from  trim  (QR)  is  assumed  to  be 
proportional  to  the  change  in  collective  blade  angle  (6^.).  This 
follows  because  blade  angle  changes  will  produce  load  changes  which  will 
result  in  torque  variations.  The  gearing  K6  is  calculated  from  flight 
test  results  so  that  QR  is  given  as  a  percentage  of  'nominal  maximum' 
torque;  i.e. 


(9  ) 


The  difference  between  the  engine  torque  QR  and  rotor  torque  QR  is  Q, 
such  that: 

Q  =  Qg  -  Qg  0  2) 

Then,  relating  torque  difference  to  rotor  speed  change  (flR): 


riR  =  (  4-)  Q  =  k5q 


(1  3) 


where  IR  is  the  rotor  moment  of  inertia  about  its  shaft  (i.e.  z-axis). 
In  addition,  the  speed  change  of  the  rotor  (a  )  is: 

n 


SR  '  /  SR  dt 


(19) 


It  should  be  noted  that  rotor  speed  (angular  velocity)  in  flight  is 
measured  relative  to  the  fuselage.  However,  any  change  in  rotor  angular 
velocity  produces  a  reaction  on  the  fuselage  which  causes  the  fuselage 
to  have  an  angular  velocity  change  relative  to  a  point  fixed  in  space 
but  in  a  direction  opposite  to  that  of  the  rotor.  Hence  the  fuselage 
angular  velocity  must  be  taken  into  account  when  considering  the 
feedback  of  rotor  RPM  to  the  governor.  However: 


ft 


F 


r 


where  r  is  the  aircraft  yaw  rate  change  from  trim 


lb) 


.  .  a 


nR  *  r 


M6) 


(  10  ) 


where  a  is  the  rotor  RPM  change  from  trim  sensed  by  the  governor 
and  S1R  and  r  are,  by  definition,  of  opposite  sign. 

Flight  results  indicate  that  for  a  collective  step  input,  the 
change  in  r  is  roughly  four  times  less  than  the  change  in  0.  This 
follows  from  the  fuselage  moment  of  inertia  about  its  z-axis  being 
approximately  four  times  greater  than  the  rotor  moment  of  inertia.  In 
the  present  exercise  r  is  neglected  as  only  pure  vertical  motion  is 
being  considered,  but  the  statement  above  should  be  noted. 
Nevertheless,  the  effect  on  overall  vehicle  response  should  be  small. 
Variable  r  is,  of  course,  required  to  be  included  when  the  yaw  channel 
and  collective  -  yaw  interaction  are  considered  (reference  5). 

2.3-  Simulation 

The  specialized  language  CSMP-10  (ARL)  (reference  12)  was  used  for 
simulation  of  the  above  equations.  However,  in  order  to  obtain  actual 
values  of  variables  rather  than  simply  changes  in  value  about  the 
(hover)  trim  flight  condition,  constants  and  integrator  initial 
conditions  need  to  be  included  in  the  equations  where  outputs  of  actual 
values  are  required  for  comparison  with  flight  test  results.  These  are 
detailed  below.  A  block  diagram  for  simulation  of  the  rotor  system  is 
shown  in  Figure  4.  Notation  used  on  this  diagram  is  standard  for  the 
CSMP-1 0 ( ARL)  language  and  corresponds  closely  with  the  notation  used 


previously. 


(11) 


2.3*1 •  Flying  controls 

The  change  in  collective  pitch  angle  (THCST)  is  generated  as  a 
function  of  time  and  converted  to  collective  stick  angle  measured 
relative  to  the  aircraft  horizontal  datum  (COLL  STK)  by  the  addition  of 
a  constant.  Thus: 

COLL  STK  -  THCST  ♦  CONSTANT  (17) 

Similarly,  actual  collective  blade  angle  (THETA  C)  is  computed  from 
blade  angle  change  (THEC)  thus: 

THETA  C  -  THEC  +  CONSTANT  08) 

In  addition,  if  an  autopilot  term  is  included,  then: 

THCAP  -  (0.86)  THCST  +  AUTOC  09) 

where  AUTOC  is  the  autopilot  3tgnal  and  THCAP  is  an  intermediate 
variable  within  the  flying  controls  simulation  which  is  generated  as 
shown  on  Figure  k. 


2.3.2.  Engine 

To  obtain  actual  engine  torque  (TORQUE)  from  torque  change  (QE)  a 
constant  equal  to  the  hover  torque  needs  to  be  Included.  Hence: 

TORQUE  -  QE  ♦  CONSTANT  (20) 


4 


(  12  ) 


2.3.3.  Rotor 

Because  Q  (the  difference  between  engine  torque  and  rotor  torque) 
is  generated  as  a  percentage  of  'nominal  maximum'  torque,  the  units  of 
rotor  angular  velocity  sensed  by  tne  governor  (OMEGA)  are  non¬ 
standard.  To  obtain  rotor  RPM  (ROTOR  RPM)  from  OMEGA,  a  gain  and 
constant  need  to  be  incorporated  thu3: 

ROTOR  RPM  =  (2600)  OMEGA  ♦  CONSTANT  (21) 


2.4.  Transfer  Function  Analysis 

Using  the  system  structure  shown  in  Figure  2,  the  transfer 
functions  and  algebraic  equations  developed  in  section  2.2  can  be 
integrated  into  the  block  diagram  illustrated  in  Figure  5.  Note  that 
because  aircraft  yaw  rate  is  not  considered  in  this  exercise,  then 

flR  ■=  a.  (22) 

Suitable  values  for  constants  are: 


K, 

- 

0.86 

t. 

- 

0.1 

k2 

- 

9.09.10- 

- 

280.0 

K, 

. 

•* 

i  > 

10.0 

K„ 

- 

K,.K„  =•  7.0  x  10 

T- 

- 

0.275 

k5 

- 

9.26. 10-* 

Is 

0.075 

K. 

- 

10.57 

If 

- 

0.26 

For 

the  purpose  of  analysis, 

,  it 

is 

convenient  to 

calculate 

the 

transfer 

functions  flR/6CST  and  Q/9C> 

the 

latter  being 

of  use 

in 

* 


computing  the  transfer  function  of  the  complete  aircraft  vertical 
dynamics  (see  section  3.9).  By  combining  blocks  and  using  the 
relationships: 


it  can  be  shown  that: 

s’CK^TjTjTj-K,  ,T,ThT6]  *Si[K,0(T,T2«-T1T5  +  TjT5)  *K  2  v  (  T  ,  T  „  +T  ,T  6  +  T  ,  T  «  )  ] 

nR  +3[K10(t,*t2  +  t5)  -  X ,  |  (Tj  +  t»  +  t6)]  +  [K10  -  Km] 

0CST  T ,T, (s+1/t , ) (s+1/t J )[s3 (t»ts )  +s2(t,+t6)  *S ( 1 +K, T s ) +K, ] 

(2*0 


On  substituting  the  constants  above: 


“r _ _  -  (o.ooeessxs1  ♦  3.5s2  -  16.0233  *  1.538) 

0CST  (s  +  0. 1 ) ( s  +  10)(s>  ♦  7 . 4285s2  *  20.785s  +  90.657) 


and  factoring  the  cubic  polynomials  gives: 


°R  _  -  (0. 006685) (s  ♦  6. 1 5) (s  -  2.5fl)(s  -  0.1) 

"  (3  +  0. 1 ) (s  ♦  t0)(s  +6.H«)(s  +  0.52  ±  J3-72) 


(  14  ) 


Similarly,  it  can  be  shown  that: 

Q _  -  (72. 1 9) ( a) (a  ♦  6.15)(a  -  2.54)(s  -  0.1) 

®CST  *  (s  +  °-1)(s  +  i°Ms  +  6.44)(s  +  0.52  ±  J3.72)  {  U 

and 

Q_  -(8.39)(s)(s  +  6.1 5) ( s  -  2.54)(s  -  0.1) 

8  '  (s  +  0.1 ) (3  +  6 . 44 )  (s  *  0.52  ±  j3 -72 )  1  0J 

The  above  transfer  functions  are  useful  for  stability  analysis. 

3.  AIRCRAFT  VERTICAL  MOTION  DYNAMICS 
3.1.  Basic  Mathematical  Model 

From  reference  7,  the  linearized  equation  of  motion  for  the 

vertical  dynamics  of  an  SH-3D  (Sea  King)  Helicopter  is  given  as: 

w -  -(0. 0001 1 )u  + (0 .000036) v  -(0.3242)w  ♦ (0 .0000057) J+ (0 .0001 8 )8-(0 . 424) 0C 

(29) 


where 


w  *  normalized  vertical  velocity  ;  change  from  trim 

u  =■  normalized  longitudinal  velocity  ;  change  from  trim 

v  -  normalized  lateral  velocity  ;  change  from  trim 

i  =  roll  rate  ;  change  from  trim 

9  =  pitch  rate  ;  change  from  trim  and 

9C  -  collective  blade  angle  ;  change  from  trim 


In  addition,  from  section  2.2.1,  the  flying  controls  equation  including 


the  autopilot  term  is: 


(  15  ) 


ec  -  <8.6)8cst  -  (10.0)6C  *  CKca)6Ca  (30) 


% 


where  -  collective  autopilot  demand  signal;  change  from  trim  and 
kCA  "  collective  autopilot  gain. 

A  block  diagram  representation  of  these  equations,  which  uses 
transfer  function  notation  and  includes  cross-couplings,  is  shown  in 
Figure  6. 

If  the  cross-couplings  are  removed  from  equation  (29),  the 
relationship  between  8C  and  w  has  the  form: 


where  z6c  is  the  gain  between  8C  and  w.  The  form  of  this  equation  is 
the  same  as  that  specified  by  Heffley  (reference  13).  The 
term  (s-z^)  is  attributed  to  the  heave  damping  mode  and  the  value 
of  zy  in  hover  is  inversely  proportional  to  the  square  root  of  the  disc 
loading.  From  Wilson  (reference  19): 

zw  ’  "  8  ^  2W/A  C32) 

where  g  =  gravitational  acceleration 

p  «  air  density 

W  «  aircraft  weight  and 

A  -  disc  area 


(  16  ) 


and  substituting  data  for  Sea  King  at  the  hover  weight  used: 


z  =■  -O.itil 

w 


This  compares  with  z  ■  -0.32142  for  the  SH-3D  given  by  Murphy  ana 

w 

Narendra  in  reference  7. 

Initially  a  computer  model  using  Murphy  and  Narendra 's  equation 
together  with  a  representation  of  the  flying  controls  was  constructed. 
In  order  to  compare  this  with  the  flight  test  results  of  reference  6,  a 
common  collective  stick  step  input  was  used.  Comparison  of  time 

histories  of  the  variables  0CST>  6C,  w  and  w  are  shown  in  Figure  7.  As 
can  be  seen,  the  quantities  w  and  w  for  the  model  bear  no  resemblance  to 
those  for  the  aircraft. 

3.2.  Model  Development 

In  order  to  overcome  some  shortcomings  of  the  basic  mode.,  the 
gain  z6c  was  increased  from  0.142k  to  3.2b  (Figure  3).  However,  there  is 
no  clear  reason  why  such  a  large  change  in  z6c  is  required. 

Nevertheless,  when  incorporated  this  gives  rise  to  reasonable  agreement 

in  w,  although  comparison  of  the  behaviour  for  w  is  poor  for  about  1 .Ss 
after  application  of  the  input. 

To  achieve  improved  agreement  between  model  and  flight  results 
for  a  'blade  angle  derivative  term'  was  added  to  the  model.  If 
cross-coupling  terms  are  omitted,  the  equation  relating  vertical 
velocity  to  blade  angle  then  becomes: 


(  17  ) 


► 


i 


> 


w  -  -(0.3242)w  -  (3.2t,)  e_  -(0.92)0. 


where  9Q  «  8r  -  0.1  0D  ,  3b ; 

Time  histories  for  the  variables  0„  ,  e^.w  ana  w  are  shown  in  Figure 
9.  It  will  oe  seen  that  the  initial  response  in  vertical 
acceleration  w  is  now  approximated,  although  behaviour  of  the  model 
after  this  initial  pulse  still  requires  improvement.  The  reason  for  the 
high  Initial  acceleration  is  not  understood,  although  a  similar 
phenomenon,  albeit  less  pronounced,  occurs  in  other  flight  results  ie.g. 
references  1  and  3).  As  rotor  speed  remains  reasonably  constant,  it  is 
possible  that  the  initial  pulse  in  w  arises  from  a  change  of  pitch.  It 
is  postulated  that  when  the  blade  angle  is  changed  rapidly  as  in  the 
step  case  used  here,  an  incidence  angle  change  exists  for  a  short  time 
until  steady  inflow  is  re-established.  During  this  transient  a 
differing  level  of  thrust  may  be  produced  which  gives  rise  to  the  change 
in  w.  However,  the  reasons  for  the  high  initial  acceleration  require 
further  investigation,  including  the  incorporation  of  improved  flight 
test  instrumentation  (e.g.  direct  blade  and  coning  angles),  so  that  a 
better  data  base  is  establisned. 

As  a  final  step  In  determining  the  structure  of  the  model,  it  was 
noted  that  flight  results  for  w  exhibited  a  periodic  variation  following 
the  initial  Impulse.  Observations  showed  that  the  magnitude  of  this 
periodicity  was  proportional  to  the  diffe"ence  between  the  rotor  and 
engine  torques,  and  that  it  was  al3o  in  phase  with  this  torque 


(  is  ) 


difference.  A  similar  periodic  variation  in  normal  acceleration  ;s 
illustrated  in  reference  3.  and  it  would  not  appear  unreasonable  trial 
acceleration  and  torque  difference  were  related.  On  tr.e  otner  hand, 
both  torque  difference  and  w  variations  may  be  due  to  blade  lag  combined 
with  pitch/lag  coupling.  Again,  improved  flight  test  information  or, 
rotor  behaviour  is  required.  However,  for  expediency  it  has  been 
decided  to  use  the  empirical  approach  here  rather  than  attempt  to 
incorporate  such  terms.  Nevertheless,  this  is  another  area  wr.ere 
further  work  is  required. 

Figure  10  shows  time  histories  of  variables  for  a  collective  step 
down  from  the  hover  when  a  torque  difference  term  is  added  to  the 
vertical  motion  equation,  while  Figure  11  is  similar  but  for  a 
collective  step  up.  In  both  cases  time  histories  for  engine  torque  ( QE ) 
and  rotor  RPM  are  included,  the  rotor  system  dynamic  modei  being  that 
described  in  section  2.  It  will  be  noted  that  the  time  histories 
for  w  are  improved  from  previous  results.  In  neither  case  have  the 

model  parameters  been  optimized  to  obtain  the  results,  as  this  would 

best  be  done  using  parameter  estimation  techniques. 

In  its  final  form,  the  vertical  equation  of  motion  for  the  model 
is : 

w  -  -(0.00011  )u  *  (0 .000036)  v  -(0.32*t2)w  *(0. 0000057)b+0. 00018)6-,  3- 2b)6c 
-(0.92)8  +(0.15)Q  18b) 

where  Q  is  defined  in  equation  (12).  In  addition,  the  constants 

associated  with  the  cross-coupling  terms  have  not  yet  been  verified 


against  flight  results. 


(  19  ) 


3.3-  Simulation 

The  simulation  block  diagram,  in  CSMP-1 0(AHL/  format,  for  vert. cal 
motion  with  the  model  in  its  finai  form  is  shown  in  Figure  t ..  ,  This 
includes  representation  of  autopilot  and  cross-coupling  terms  together 
with  means  for  obtaining  actual  values  of  variables  as  describe.:  iti 
section  2.3.  WKN  (vertical  velocity  in  Knots.'  is  istair.eo  from 

W( vertical  velocity  in  ft/s)  and  altitude  HT  is  obtained  by  integrating 
W.  The  model  time  histories  shown  in  Figures  10  and  '1  use  this 
simulation  combined  with  that  shown  in  Figure  reasonable  comparisons 
of  the  time  histories  of  variables  have  also  been  obtained  for  pulse  and 
doublet  collective  stick  inputs. 

3.*i.  Transfer  Function  Analysis 

When  the  autopilot  input  and  cross-coupling  terms  are  omitted  from 
the  final  version  of  the  model,  the  transfer  function  block  diagram  li¬ 
as  shown  in  Figure  13.  The  transfer  function  for  trie  flying  controls 
(  is  33  specified  in  equation  (2)  but  the  vertical  motion 

dynamics  ( w/ Wq )  now  comprise  the  following: 

i.  The  blade  angle  derivative  term 

6CP  =  9.2s  ,  7 

ec  s  +  10 

which  can  be  obtained  from  equations  (3'*j  and  (3b). 


ii.  The  proportional  term 


r 


i 
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» 


\ 


°c 

which  is  simply  the  increased  gain  z6c  and 

iii.  The  torque  difference  term 

9CQ  1  ■  26(  3  )  ( s  ♦  6 ■  1  5)  (3  -  2  ■  5*0  ( s  -  0.1) 

0C  (s  ♦  0. 1 )  (s  -  6.44)(s  *  0.52  ±J3.72)  U 

which  i 3  derived  from  equation  (28)  (the  Q/9(.gT  transfer  function 
multiplied  by  the  appropriate  gain). 


The  vertical  velocity  can  then  be  computed  based  on  equation  (31) 
thus: 


w  _  _ 1 

8C'  "  3  +  0.3242 


(40) 


where  9C' 


CD 


CP 


*  8 


CQ 


(41) 


It  should  be  observed  that  both  w  and  h  are  measured  positive  downwards 
in  standard  notation  and  that  height  h  is  obtained  from  w  thus: 


h 


w 


(42) 


If  all  terms  are  included,  the  overall  transfer  function  is: 


h  -  117.9(s  +  3. 2) (s  *  0.36  ±J 3 . 2 ) 

0CST  *  ( 3 )  ( s  <■  0.  3242)  ( 3  *  1  0)  (s  *  ToHs  +  0.52  ±  j 3. 72) 


k 


i 


(43) 
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However,  if  some  slmpl  if  teat  ion  is  required  for  analysis  purposes  where 
higher  order  terms  are  not  Important,  the  torque  difference  term  can  De 
omitted  to  give: 

h  -  1  07 ( s  ♦  2.61) _ 

ecST  ’  (s)(s  +  1 0)  (s  *  1  0 )  ( s  +0.  32142)  v  1 

or,  with  the  derivative  term  also  omitted: 

.  _ zJldl _  ,„5) 

0CST  (s)(s  ♦  1 0)  (s  ♦  0.3242)  ' 

3.5.  Stabilization 

This  section  serves  to  illustrate  how  the  transfer  functions 
generated  in  section  3.4  can  be  used  when  height  stabilization  systems 
are  incorporated  in  the  helicopter.  The  Sea  King  Mk  50  helicopter  is 
fitted  with  two  forms  of  height  stabilization,  namely  Radio  Altitude 
Hold  and  Barometric  Height  Hold.  The  former  provides  accurate  height 
stabilization  up  to  500  ft  and  also  enables  automatic  transition  and 
hover  manoeuvres  to  be  performed.  Barometric  height  hold  provides 
stabilization  in  height  up  to  ceiling  altitude.  Full  details  of  both 
systems  are  given  in  references  15  -  IT. 

3.5.1.  Barometric  height  hold 

In  order  to  configure  the  system  in  a  manner  suitable  for  stability 
analysis,  it  is  necessary  to  make  a  number  of  assumptions  and 
simplifications  to  the  actual  aircraft  flight  control  system.  These 


include: 


(  22  ) 


i.  Ignoring  ’open-loop  spring  operation’  (see  reference  15)  as  this 
serves  merely  to  re-datum  the  system  when  relatively  large 
autostabilizer  signals  are  generated  and  is  non-linear  in  its 
operation. 

ii.  Ignoring  the  (non-linear)  authority  limits;  i.e.  only  small 
perturbations  about  a  given  flight  condition  are  considered. 


When  these  simplifications  are  made,  the  control  law  governing 
Barometric  Height  Hold  is  (from  reference  15): 


,  0.0622 

6CA  '  '  s  -  2.22  J(h  "  "ref1 


where  hpef.  is  the  reference  altitude  and  the  smoothing  term  (lag)  is  a 
gust  filter.  A  block  diagram  incorporating  this  law  and  the  aircraft 
vertical  motion  dynamics  is  3hown  in  Figure  1^.  This  forms  the  basis  of 
analysis  to  determine  barometric  height  hold  system  stability. 


3.5.2.  Radio  altitude  hold 

As  with  Barometric  Height  Hold,  it  is  necessary  to  configure  the 
system  in  a  manner  suitable  for  stability  analysis.  The  Radio  Altitude 
Hold  circuit  used  in  the  actual  aircraft  flight  control  system  includes 
an  inertial  height  smoother  and  is  quite  complex.  Both  the  system  and 
its  mathematical  model  are  described  in  reference  IT.  A  simplified 
version  of  this  i3  presented  in  reference  15,  although  even  Lhi~ 
contains  non-linear  elements.  Consequently  simplifies ~ 
assumptions  similar  to  those  used  in  section  3*5-1  need  to  be 
obtain  a  representat ion  suitable  for  stability  analysis. 


When  simplifications  are  incorporated  in  the  model  given  in 


reference  15,  the  control  law  governing  Radio  Altitude  Hold  becomes: 


„fS2  *  0.22s  +  0.036, 
9CA  =  Ch  href)|0‘371  sfs“  0.22] - J 


,  ,1  0.062  :  . 
h  o.'as1  *  °'"ow  “ 


and 


HS 


s  +  T 


HS 


A  block  diagram  representing  this  law  is  shown  in  Figure  "5  ana  again 
this  forms  the  basis  for  system  stability  analysis.  It  should  be  noted 
that  an  inertial  height  smoother  modifies  the  height  signal  n  via, 
basically,  a  simple  lag  ci-cult,  to  give  h'.  Provision  for  tnis  :s 
indicated  on  Figure  15,  where  KHS  and  are  both  estimated  to  te  1.3. 


4.  CONCLUDING  REMARKS 

A  mathematical  model  representing  the  dynamic  flight  behaviour  of  a 
hovering  helicopter  has  been  constructed.  Only  vertical  motion  at  ut 
the  hover  point  has  been  considered.  The  model  developed  comprises  two 
parts,  the  rotor  system  dynamics  and  the  aircraft  vertical  motion 
dynamics.  To  obtain  an  accurate  representat ion  of  vertical  fiignt 
behaviour,  it  has  been  found  necessary  to  incorporate  a  rotor  system 
term  in  the  flight  dynamic  model  for  vertical  motion. 

The  rotor  system  model  consists  of  a  number  of  components  which  are 
linked  together.  A  dynamic  representation  of  each  element  has  been  made 
and  parameters  adjusted  so  that  the  dynamics  of  the  model  are  similar  to 
those  of  the  physical  system.  Simulation  of  the  model  is  discussed  and 
a  transfer  function  analysis  carried  out. 


(  2*i  ) 


The  aircraft  vertical  motion  dynamic  model  is  based  on  previously 
published  work  which  includes  cross-coupling  terms.  However,  in  order 
to  obtain  an  accurate  representation  of  aircraft  vertical  motion,  a 
number  of  additions  to  the  basic  model  have  been  found  necessary.  In 
conjunction  with  adjusting  the  gain  relating  vertical  velocity  change  to 
collective  blade  angle  change,  extra  terms  have  been  added.  These  are  a 
blade  angle  derivative  term  and  a  torque  difference  term  which  are 
required  so  that  the  model  behaviour  can  be  adequately  fitted  to  the 
flight  test  results  for  similar  manoeuvres.  Whilst  some  bases  for  these 
additional  terms  have  been  postulated,  further  work  is  required  to 
establish  a  better  understanding  of  flight  behaviour.  This  includes 
instrumentation  of  additional  rotor  variables  during  flight  test. 

After  these  modifications  are  made,  the  flight  dynamic 
characteristics  of  the  model  are  shown  to  be  similar  to  those  of  the 
flight  test  results.  However,  no  attempt  has  been  made  to  optimize 
parameters  as  this  would  best  be  done  using  identification  methods.  The 
linear  model  developed  is  ideally  suited  to  this  approach.  Again, 
simulation  and  transfer  function  analyses  have  been  conducted;  the 
latter  is  of  use  in  determining  both  the  stability  of  the  bare  aircraft 
and  system  stability  when  height  hold  loops  are  added. 

The  technique  used  to  develop  the  model  structure  has  been 
demonstrated  to  give  a  realistic  representation  of  vertical  flight 
behaviour  for  small  perturbations  from  the  hover.  The  derived  transfer 
functions  should  be  of  use  for  handling  quality  and  flight  control 
system  studies  and  the  linearised  model  should  be  useful  for  gust 
response  analyses.  It  is  proposed  to  use  the  technique  to  observe 
flight  behaviour  in  other  axes,  both  in  hovering  and  forward  flight,  so 
that  a  greater  understanding  of  helicopter  dynamics  can  be  obtained. 


(  25  ) 


REFERENCES 

1.  Wilcock,  T.  and  Thorpe,  A.C.  :  Flight  Simulation  of  a  Wessex 

Helicopter  -  A  Validation  Exercise.  RAE  Tech.  Report  73096  (1973). 

2.  Wilcock,  T.  :  A  Piloted  Flight  Simulation  of  the  Westland  lynz.  RAE 
Tech.  Report  79099  (1974). 

3.  Padfield,  G.D.  :  A  Theoretical  Model  of  Helicopter  Flight  Mechanics 
for  Application  to  Piloted  Simulation.  RAE  Tech.  Report  TR  81048 
(1981). 

4.  Guy,  C.R.,  Williams,  M.J.  and  Gilbert,  N.E.:  A  Mathematical  Model 
of  the  Sea  King  Mk  50  Helicopter  in  the  ASW  Role.  ARL  Aero  Report 
156  (1981). 

5.  Guy,  C.R.:  Hovering  Helicopter  Flight  Dynamics  -  A  Study  of  Yawing 
Motion.  ARL  Aero.  Tech.  Memorandum  to  be  published. 

6.  Guy,  C.R.  and  Williams,  M.J.:  Sea  King  Helicopter  Flight  Trials. 
ARL  Aero.  Note  415  (1983). 

7.  Murphy,  R.D.  and  Nardendra,  K.S.:  Design  of  Helicopter 

Stabilization  Systems  Using  Optimal  Control  Theory.  J.  Aircraft, 
Vol.  6,  No.  2  (1969). 


(  26  ) 


8.  Packer,  T.J.:  Wessex  Helicopter/Sonar  Dynamics  Study:  The 

Mathematical  Model  of  the  Helicopter  Ae-odynamics  and  Kinematics. 
WRE  Tech.  Memorandum  SAD203  0  969). 

9.  Guy,  C.R.:  Sea  King  Mk.50  Helicopter  Flight  Control  System:  A 

Mathematical  Model  of  the  Flying  Controls.  ARL  Aero.  Note  388 
0979). 

10.  Blakelock,  J.H.:  Automatic  Control  of  Aircraft  and  Missiles. 

Wiley  0  965). 

11.  :  Course  Notes  on  Sea  King  Helicopter  Engine  Control  System. 

12.  Gilbert,  N.E.  and  Nankivell,  P.G.:  The  Simulation  Language  CSMP-10 
(ARL).  ARL  Aero  Note  362  0976). 

13.  Heffley,  R.K. :  A  Compilation  and  Analysis  of  Helicopter  Handling 
Qualities  Data.  Volume  Two:  Data  Analysis.  NASA  Contractor  Report 
31  45  0979). 

19.  Wilson,  E.B.:  Theory  of  an  Airplane  Encountering  Gusts,  II,  NACA 
Report  No.  21  0916). 

15.  Guy,  C.R.:  Sea  King  Mk.  50  Helicopter  Sonar  Dynamics  Study.  A 
Simplified  Control  Systems  Mathematical  Model.  ARL  Aero.  Report 


152  (1979). 


(  27  ) 


16.  Guy,  C.R.:  Sea  King  Mk.  50  Helicopter  Flight  Control  Systems.  A 
Mathematical  Model  of  the  AFCS  (Autostabilizer/Autopilot  Mode). 
ARL  Aero.  Note  387  (1979). 

17.  Guy,  C.R.:  Sea  King  Mk  50  Helicopter  Flight  Control  System.  A 
Mathematical  Model  of  the  AFCS  (ASW  Mode).  ARL  Aero.  Note  393 


(1979). 


Collective  stick 
angle  (deg) 


0 

10 


Collective  blade 
angle  (75%  radius) 
(deg) 


Normal  acceln. 
(ft/-,*) 


Torque  (1) 
(percent) 


Rotor  rpm 
(percent) 


FIG.  1.  COLLECTIVE  DOUBLET  INPUT  -  80 
(FLIGHT  TEST  RESULTS) 
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COMPARISON  OF  MODEL  AND  FLIGHT  TEST  RESULTS  (BASIC  MODEL) 
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COMPARISON  OF  MODEL  AND  FLIGHT  TEST  RESULTS  (INCREASED  GAIN 
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COMPARISON  OF  MODEL  AND  FLIGHT  TEST  RESULTS  (DERIVATIVE  TERM 
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FIG.  10.  COMPARISON  OF  MODEL  AND  FLIGHT  TEST  RESULTS  (TORQUE  TERM  ADDED) 
-  COLLECTIVE  STEP  DOWN 


AIRCRAFT  VERTICAL  MOTION  TRANSFER  FUNCTION 


FIG.  14.  BAROMETRIC  HEIGHT  HOLD  BLOCK  DIAGRAM 
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